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— . Abstract 

qh| The classical electromagnetic modes outside a long, straight, su- 

perconducting cosmic string are calculated, assuming the string to 
bX)' be surrounded by a superconducting cylindric surface of radius R. 

Thereafter, by use of a Bogoliubov-type argument, the electromag- 
netic energy W produced per unit length in the lowest two modes is 



^ ' calculated when the string is formed "suddenly". The essential new 

^ element in the present analysis as compared with prior work of Parker 

[Phys. Rev. Lett. 59, 1369 (1987)] and Brevik and Toverud [Phys. 
Rev. D 51, 691 (1995)], is that the radius a of the string is assumed 
finite, thus necessitating Neumann functions to be included in the 
fundamental modes. We find that the theory is changed significantly: 
Wis now strongly concentrated in the lowest mode (m, s) = (0,1), 
whereas the proportionality W oc (Gfi/t)'^ that is characteristic for 
zero-width strings is found in the next mode (1,1). Here G is the 
gravitational constant, the string mass per unit length, and t the 
GUT time. 
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1 Introduction 



Cosmic strings are remnants of an ultrahigh-temperature phase transition 
in the early universe, allowed in certain types of field theories [|l|. They 
are explained in general terms of spontaneous symmetry breaking and have 
analogies to vortex lines in superfluid helium and in superconductors. If they 
exist, they may have seeded galaxy formation and may explain the observed 
fluctuations in the cosmic background radiation. 

An interesting aspect of cosmic string theory is the estimate of the amount 
of radiation produced when the string is formed " suddenly" , at some instant 
t = to, where to — 10~^^s is the characteristic time for the grand unified field 
theory (GUT). One can express the number of massless particles created by 
the formation of the string in terms of Bogoliubov coefficients relating the 
initial metric to the (static) metric after the string is formed. This kind of 
approach was pioneered by Parker 0, for the case of a massless scalar field. 
As the radius a of a string is believed to be very small (~ 10"^^ cm for 
symmetry breaking at the GUT scale) the string was in the formalism put 
equal to zero, thus making it possible to avoid the terms in the solution of 
the field equations that diverge at the centre (r = 0). If factors of unity are 
disregarded, the analysis of Parker shows that the energy produced per unit 
length in the lowest non vanishing radiation mode is of order W ~ (G/x/t)^, 
where fi ~ 10^^ g/cm is the string mass per unit length (GUT scale). 

An analogous calculation to that of Parker was made in Ref. in the 
electromagnetic case, still under the assumption that the string radius a is 
equal to zero. This calculation is pertinent to the particular subclass of su- 
perconducting strings , believed to be existent in the early Universe under 
certain conditions. We shall not enter into a discussion of the superconduct- 
ing string's physical properties here. Rather, we take a pragmatic attitude 
and adopt a physical model in which the string can be taken to obey perfect 
conducting boundary conditions at r = a. This kind of model is motivated 
by the properties known from electromagnetic theory for a superconductor 
0: The magnetic field must be everywhere tangential to the surface, and 
the electric field must be everywhere normal to the surface. Thus we have 
precisely the same boundary conditions as if the surface were perfectly con- 
ducting. 

In 0] we found essentially the same result as Parker did in the scalar case: 
For the lowest non vanishing mode, the energy W was found to be of order 
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The purpose of the present paper is to calculate, again with use of the 
qnantTim mechanical sudden approximation, the amount of electromagnetic 
energy produced if the string has a finite radius a. The limit a — > is more 
delicate than one might anticipate at first: The Neumann functions present 
in the formalism for finite a blow up in this limit, and one cannot in advance 
know what the physical implications of this are, whether the emitted energy 
W diverges, stays finite, or goes to zero. A concrete calculation is needed. 

In the next section we give the exterior Gott-Hiscock metric, which is a 
central element in the subsequent analysis. Thereafter, in Sec. 3 we calculate 
the three classical modes in the cavity, and the energy associated with the 
TM and TE modes. In the quantum mechanical theory in Sees. 4 and 5 
we consider the TM mode only, and calculate by means of the Bogoliubov 
transformation the electromagnetic energy Wms produced in the mode (m, s), 
for a GUT string. The energy production in the lowest mode (0, 1) is found 
to be large though finite: Wqi ~ 6 x 10^^ erg/cm. It is not until the next 
mode (1,1), where the energy is much smaller, that W is found to obey the 
proportionality mentioned above: Wn oc {Gn/ty. The inclusion of the finite 
string width thus changes the theory in a significant way. 

We do not consider the limit a explicitly in this paper. We calculate 
however the emitted energy for a relative string thickness in the interval 
10~^ < a/R < 10~^. The result is that over this range the energy in the 
(0,1) mode is still dominant. 

We adopt in the following Heaviside-Lorentz units, and put h — c— 1. 
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2 The string metric 



In this section we collect some information about the metric around an in- 
finitely long, straight string of finite radius a. Let the string's centerline 
coincide with the z-axis. We employ cylindrical coordinates, numerating 
them as = {t, r, 6, z). 

We first note that within the conventional approach, when neglecting the 
string's width and using the linearized Einstein's equations, one has for the 
exterior metric 

dsl^^ = -dt^ + dr^ + (1 - SGfxyde^ + dz^, (1) 

where G is the gravitational constant and fi, as mentioned, is the string mass 
per unit length. The solution (1) was found by Vilenkin, within perturbation 
theory 0. The deviation from the Minkowski metric is small, Gfi ~ 10~^. 

Then let us take into account the finite radius a. We assume the mass 
density p(r) to be uniformly distributed across the string, 

p(r) = I ' ^ "' (2) 

'^^ ^ [ for r > a. ^ ^ 

Then, with the abbreviation 

^={8nGpoY/\ (3) 
we can write the interior metric as 

dsl, = -dt" + dr^ + (^^^^\e^ + dz\ (4) 

where — oo < t < oo, < r < 9m/i^ where 9m is a constant, < 9 < 2n, 
and — oo < z < oc. 

To find the metric on the outside, we have to take into account the con- 
ditions that the metric itself, as well as its first derivative with respect to r, 
have to be continuous across the surface of the string (these are the Israel 
matching conditions 0). We get 

dsl^, = -df + dr^ + (1 - 4Gfifr^d9^ + dz\ (5) 

This is the exact metric, valid for r > a, and will be used in the following. 
The constant 9m introduced above is related to the mass per unit length by 
the relation 

/i= ^(1-COS^Af). (6) 
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Strictly speaking, the metric (5) applies only for Gjj, < |. But, as noted, 
this condition is amply satisfied for cosmic strings. 

In the case of a weak field, <^ 1 which implies <^ 1, the exact 
solution (5) is seen to reduce to Vilenkin's singular line metric (1). Then Eq. 
(6) yields 8G/x = 9\j which, together with = na^pQ, yields 9m = i^ct. The 
upper limit 6m/ for r in the interior metric (4) thus becomes equal to a. 

The exact solutions (4) and (5) were found by Gott and by Hiscock 
. Cf . also the papers of Helliwell and Konkowski |T0| , Jensen , and the 
book of Peebles pi. 



3 The classical electromagnetic modes 

The superconducting string of radius a, as mentioned, is assumed to be 
formed suddenly at the instant t = to, of the order of the GUT time. We 
assume that the electromagnetic field, at all times, is constrained to lie within 
a large cylinder of length L and radius R, whose centerline is coinciding with 
that of the string (the 2;-axis). Here R and L are taken to be of the same order 
of magnitude as the mean spacing between strings at the time of formation. 
This length can be no larger than the horizon size, which in a radiation 
dominated universe is Rp = 2t. We will use R = t. For the electromagnetic 
fields we assume perfect conductor boundary conditions, on the surface r = R 
as well as on r = a. 

Prior to the formation of the string we have the Minkowski metric 

ds^ = -dt"^ + dr^ + r^de^ + dz\ (7) 

where t < to, < r < i?, < ^ < 27r, -^L < z < \L. After the 
formation, the metric is the Gott-Hiscock metric (5) in the exterior region. 
It is convenient to introduce the symbol 

f5={l-AGiiy^. (8) 

The field modes calculated for t < to are related to the modes for t > to by a 
Bogoliubov transformation. The change in the metric has the effect that the 
vacuum states associated with the two sets of modes differ. As a consequence 
there will be a net particle production in the t > to modes, relative to the 
t < to vacuum. This is precisely the effect that is calculable in terms of the 



Bogoliubov coefficients of the transformation |T3|. Before embarking on this 



dynamical calculation, we consider first the stationary classical modes. 



5 



3.1 Modes for t > to 



Taking the time factor to be exp{—iut) in a complex classical representation 
for the fields, we get Maxwell's equations in the annular region between r = a 
and r = R: 

V ■ E = 0, V ■ H = 0, (9) 



V X E = icoB., V X H = -icoB. (10) 

Expressed on component form, these equations are the same as Eqs. (5) - 
(8) in P, and will not be repeated here. From these equations we can derive 
a second order equation for Ez- We write this field component in the form 
Ez{r, 6, z) = Ez{r) expi{m6 + kz), where m is an integer, and k is the axial 
wave number. We assume that k is continuous for — oo < A; < oo. The 
equation for E^ becomes 

Ez{r) = 0, (11) 

so that the solutions are linear combinations of ordinary Bessel functions 
J|/3m|(gr) and Neumann functions N\ism\{<ir), with = uj^ — k"^. 

TM mode: 

According to the boundary conditions, E^ vanishes for r = a,R. From 
this we derive the condition 

J|,|(gi?)iV|,|(ga) - iV|,|(gi?) J|,|(ga) = 0, (12) 

where we have defined u as 

u = (3m. (13) 

For a fixed value of m, Eq. (12) gives an infinite sequence of roots q^s, with 
s = 1, 2, 3, ... As is real, we can now write 

^.sk = {ql + k'f". (14) 



(P 1 d 
dr'^ r dr 
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We let A^sk be an arbitrary mode-dependent constant, define Ci^\ as 

C\u\iqusr) = J\u\{qvsr) - ^M^^Ar|^|(g^,r), (15) 

^^\u\{qusa) 

and use the real representation for the azimuthal 6 direction. The field 
components can then be written as 



ik , 



/ cos\ 
\sm J 



Er = A.sk—Cl^iiqusr) { )m^e**, (16) 

qus 

Eg = TA.sk'^C\,\{qusr) f '''' V^^'^ (1^) 
ql,r \cosJ 

( COS \ 

= ^.sfcC|,i(g,,r) . m^e**, (18) 



n.. = ±A...^CH(g..r) ('^]m9e^'', (19) 

qLr \cosj 

He = A^J^q^^iq^sr) Mm^^^*' (2°) 

where ^ = kz — uj^skt- Prime means differentiation with respect to the 
whole argument. Note that the m = mode is independent of the azimuthal 
direction 6. 

Comparison with the analogous expressions given in [Q] shows that the 
difference between the two cases lies in the radial functions. Explicitly, 
J\,^\ C\u\, which in turn implies the presence of the Neumann function. 
The physical reason for this is, of course, the finiteness of the string's width. 

TE mode: 

This mode is constructed by starting from the governing equation for the 
component Hz- The boundary conditions now lead to the condition 

J[,\{q'R)Nl^\iq'a) - N[^^{q' R)J[,^{q' a) = 0, (21) 

which has roots qlg,s = 1,2,3... The complete expression for the TE field 
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components are 



i^. = ^s.^C[,|(g:,r)Mm^e^*', (22) 
Qi/s \smj 

He = TA.s^j^^CMsr) ('1")-^^^"', (23) 

= A,,,q,\{ql^r) Mm^e^*', (24) 
ysm J 

Er = tKsU^^C^^HsT) f '''' V^^^*'' (25) 
[ql.Yr \cosJ 

Ee = -i...^C[,|(gLr) iZV^'''" ^ (26) 
where the A^sk are a new set of constants, tol^^ = [(g^^^)^ + k'^Y^'^, and $' = 

TEM mode: 

In this mode both components and are equal to zero. The field 
E(r, 6) is the solution of a two-dimensional electrostatic problem, just as in 
the Minkowskian case |14|. With E(r, 6) = — V0 inserted into V ■ E = we 
get the generalized Poisson equation 

The boundary conditions Er{a, 9) = Er{R, 6) = imply azimuthal indepen- 
dence of 0, and we get the solution of Eq. (27) in the form (p = A — Blnr 
with A, B as constants. From Maxwell's equations we get the mode solutions 

E = :^e**f, H(r) = ^e'^'O. (28) 

It turns out that the TEM mode is special and poses problems when devel- 
oping a manageable quantum formalism for the field modes. This mode will 
therefore not be considered in the quantum theory below. 
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3.2 Modes for t < to 

The initial conditions are those of an empty cyhnder of radius R in Minkowski 
space. We do not give the complete solution here. As in the case considered 
in , the only components needed for the calculation are and Hz for the 
TM and TE modes, respectively. For t < to these are 

(COS \ 
. m^e'*, TMmode (29) 
sm J 

where p^s, s = 1, 2, 3, ... and p'^^, s = 1, 2, 3, ... are roots of 

J|^l {pR) = and J|'„| {p'R) = 0. (31) 

The TEM mode does not exist for t < Iq. This is just similar to the 
Minkowskian case. 



3.3 Energy considerations 

We note first the following point, without going into detailed considerations: 
Using Maxwell's equations, we can find out whether there is an axial electric 
current flowing in the string, in the TE and TM modes. Similarly we can 
find out whether there is a charge Q/L per unit length. Substituting the 
appropriate field components into the expressions for and Q/L we find 
that it is only the mode m = that gives an answer different from zero. 
A net current, and an electric charge, are thus associated with azimuthal 
symmetry only. 

We next concentrate on the total electromagnetic energy per unit length, 
W, outside the string: 

W = ] Td^l' + \H\'')dA, (32) 

4 Ja 

where dA = (3~^rdrd9 is the cross-sectional area element of the guide. We 
get 

^\VtEz\'^dA, TMmode 
W =\ " (33) 
[li^S\"^tHz\''dA, TEmode 
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where Vt is the transverse gradient which for the TM mode is defined by 



VtE,\'dA = - 



dr 



dd 



rdrdO. 



(34) 



Green's first identity together with the wave equation and boundary condi- 
tions for Ez gives 



j \VtE,\^dA = qlJ \E, 



(35) 



which, together with an analogous result for the TE mode, leads to the 
convenient result that the field energy is expressible entirely in terms of the 
Ez and Hz components: 



TM mode, 



W 



(36) 



/ \Hz\'^dA, TE mode. 



Explicit calculations can now be carried out using the boundary conditions 
together with the orthogonality properties of the Bessel functions |15[. We 
get for the TM mode 



W 



mis 

and for the TE mode 



'''''' ^ Aysk\^[R^C'\^\+l{(lvsR) — '^^(^|t|+l(Q'i/sO) 



(37) 



W 



A' P 



(38) 



4 Quantum mechanical modes 

As in Ref. [0], when developing the quantum formalism, we will concentrate 
on the TM mode only. We have seen that the field energy can be expressed 
solely in terms of the component Ez- We can therefore essentially make use 
of the formalism of scalar field theory. We expand Ez as a field operator: 

— I]KsfcMi.sfe(a;) + alskKsk{x)i (39) 

-OO ZTT , „ 
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where Ui^gk and al^i. satisfy the commutation relations 

Ksfc, al's'k'] = '^T^5{k - k')5^yi5ss'-, (40) 

the other commutators vanishing. Here we use as a general index, meaning 
that 

m, for t < to ^^^^ 
/3m, for t > to- 

The mode functions u,^sk are normalized according to 

■J{\E,\^)dV={a^a)u;,sk, (42) 



usk 



where dV is the volume element. This ensures that the quantum mechanical 

expression equivalent to (36) becomes equal to the expectation value of the 
number operator times the photon energy of the mode. Thus, the normalized 
mode functions for t > to are 



13 (Lini, 



xC|^^l(g^^.r)(^2'jm^e^^ (^3) 

where 

appears as a consequence of the real representation for the 6 direction. The 
corresponding functions for t < to are 



Umsk{x) = e^J PI J ^""Z ^J\m\{Pmsr)\^]m9e'^. 

]j nUJmsk R\J\m\+l{PmsR)\ J ^^^^ 

We define the Klein-Gordon scalar product as in scalar field theory: 

(01, 02) = -4 / ['^i^*'^2 - idtcPMdV, (46) 
q J 
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where — Ylivs ots- ^^^^ mode functions are then orthonormal: 

{u^sk-, K's'k') = 27r5(/c - k')6^i,'6ss', (47) 

the other products vanishing. 

For later calculations it will be important to specify the continuity condi- 
tions for the field at t = to- First, we require that the field itself is continuous: 

E,(x)\,-^E,(x)\,^. (48) 

Next, we can use the conservation of the Klein-Gordon product to find the 
continuity condition for dfEz'. 

-J / [Ez\,^{dtE:)\,^ - {dtEz)\,^E:\,^]rdrdedz = 

_. (49) 
-4 I [E^,-{dtEl)\^- - {dtE,)\^-E*X-]rdrdedz. 

Substituting Eq. (48) into this equation we get 

-,2 



d,Ez{x\- = ^^d^Ez{x\^^, (50) 



where = Em.PL- 



5 Energy production calculated from the Bo- 
goliubov coefficients 

5.1 General formalism 

The two sets of mode functions give two different expansions for the field: 
For t < to, Ez is given by Eq. (39) with p = m; ior t > to, E^ is given by 
the same expression with v — /3m. For each of these expressions there is a 
unique vacuum state: 

amsk\0)nisk = 0, for t < to, (51) 



apmsk\0)pmsk = 0, for t > to. (52) 
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As both sets of modes are complete, Ujjmsk can be expanded in terms of 

7r~ J2 Yi{l3msk\m's'k')um's'k'{x) 

27r (53) 

+ 6{(3msk\m's'k')u*^,^,f^,{x% 

7 and 5 being the Bogohubov coefficients [|l3l. Equivalently, for the operators 
we have the expansions 

TT- Yl [l{l3msk\m's'k')am's'k' 

277 (54) 

+ 5* {Idms^m' s' k')a}^, 

This imphes that the average number of particles produced in the mode 
m, s, k per unit interval in the k space is 

^ = r ^ E mnisk\m's'k')\\ (55) 
ak J-oo 2tt , I 

m' ,s' 

The coefficient 5{(3msk\m' s'k') can be found using the continuity condition 
for the field together with the orthonormality of the scalar product: 



(56) 



5{(3msk\m Sk) = / [Umskix)dtUf}rn's'k'{x) 

-(3—dtUmsk{x)upm's'k'{x)]rdrd6dz. 
Substituting Eqs. (43) and (45) into Eq. (56) we have 

6{Pmsk\m's'k') = —= ^^^^ ^ ^ ^ ( k — k') 

V P Qpm's'k' 



X 



I 2 I 

13m' s'k' nQ I ^msk 



^msk y ^Pm's'W _ 



(57) 



where 



J ^ _ la J\m\{Pmsr)C\i3m\{qi3ms'r)rdr 

\RJ\m\+l{PmsR)\\R'^Cfp„il+l{qf3ms'R) " «^(^|^/3m|+l (^/3ms'") I ^ (58) 

Our main task now is to calculate these expressions with satisfactory accu- 
racy. 
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5.2 Evaluation 



As (3 is close to unity, we can put (3 = 1 everywhere except in the difference 
between the square roots. We can also approximate the roots q^ms of the TM 
mode equation (12): As a/i? ~ 10~^ and the Neumann function diverges for 
small arguments, we can neglect the second term in the equation and so 
obtain J\/3m\{(li3mR) — 0, which means that g/j^s — Pms- Using this, together 
with the general properties of the Bessel function, we find that Eq. (58) can 
be approximated by 

Ls' ^ ^Sss'- (59) 

More detailed considerations on the evaluation of Iss' sue given in Appendix 
A. Equation (57) can now be written 



5{(3msk\m' s' k') = tt 



^msk 



^msk y'^/3m's'fc'J (60) 



This leads to the following expression for the electromagnetic energy pro- 
duction in the mode m, s, k per unit length, and per unit k interval: 

dNmsk _ 1 {^I3msk ~ (3^ msk) x 

dk L dk A (^/smsk 

As in and [0 we simplify the calculations by considering only the wave 
number interval centered around k = 0. Since Uf^msk = ilpms + k'^Y^'^ and 
t^msfc = (Pms + fc^)^^^ it follows that ujf^rnso = Qprns and Umso = Pms- Therefore, 
the energy production can be found from the roots of the Bessel function: 

dWmsk I 1 (qtSms - PPmsf 

—7, — U=o - 7 ■ (t»^J 

dk 4 qpms 

In order to make a reasonably accurate estimate of the energy production 
one has to be careful: As |1 — /9| ^ 4 x 10~^ it follows that insertion of crude 
approximations for the roots may easily lead to substantial errors. The best 
way to proceed, is to check the roots on a computer. In Table I we give a 
short compilation of roots xpms = Rqpms and Xms = Rpms calculated by use 
of Maple V. 

It is clear that in the lowest mode, (m, s) = (0, 1), the difference between 
the non gravitational and the gravitational cases is most significant. Let 
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us compare the energy produced in this mode with that of the next mode, 
(m,s) = (l,l): 



msk 



dk 



^ X 10-2 - I X 10-3, for (0, 1) 



k=o 1 1(1-/5)2- 2 ^ 10-11, for (1,1). ^^^^ 



The tendency to produce more energy in the lowest mode is very strong: The 
value for (0, 1) is 10^ times larger than the value for (1, 1). 

To get the total energy, we multiply the above equation with the width 
Ak of wave numbers contributing to the energy production. As the string 
actually forms during a finite time interval At, the production of photons 
with wave numbers much larger than 1/At will be suppressed. We take 
At = R = t. For the total energy per unit length we then obtain 



1 dWmsk 





t dk 



I X 10-3, for (0, 1) 
fc=o I^xl0-i\ for (1,1). 



(64) 



As mentioned above, we take t to be the GUT time, t ~ 10 §. For the 
(0, 1) mode we then have, in ordinary units, 

H/q^ ~ X 10-3 ~ 6 X 10^9^, (65) 
c't^ cm 

and for the (1, 1) mode, 

TV,, ~ 4^ X 10-^1 ^ 6 X 10^1^. (66) 

c^t2 cm 

Characteristic for the finite-width string considered here is that the dominant 
mode is the one corresponding to azimuthal symmetry, m = 0. For the zero- 
width string considered in 0, we got instead the emitted energy in this 
mode to be zero. Inclusion of the finite width thus changes the theory in an 
essential way. The lowest non vanishing mode found in was (m, s) = (1, 1). 
From Eqs. (63) and (64) we see that, in the present case, the (1, 1) mode 
corresponds to Wn being proportional to {G^/ty. This agrees qualitatively 
with the behavior found in ||^ for the same mode. Quantitatively, Eq. (66) 
is of the same order of magnitude as Eq. (57) in M. 
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5.3 On the a/R- dependence 



In view of the sensitivity that we have seen in the calculation of the energy 
production, we may ask to what extent the results are dependent on the 
values chosen for the input parameters. 

Recall that our input parameters are the string mass /i per unit length, 
the string radius a, the string length L, and the radius R of the waveguide. 
In terms of these parameters the energy production per unit length can be 
written as (cf. Eq. (61)): 

ALR Xpms 

Here it is to be noted that the roots xp^s depend on the ratio a/R. Table II 
gives the roots for various values oi a/R (they were calculated using Maple 
V). Note also that the roots Xms are independent of a and i?, and are as 
given in Table II. Prom these roots we have calculated the a/R- dependent 
results for energy production as shown in Table III. 

It is seen that with increasing values of the string width, the amount of 
produced energy increases. This is in accordance with expectations. The 
increase is however relatively moderate as long as the string stays "thin" 
relative to the radius of the waveguide: In the range 10~^ < a/R < 10~^ the 
energy increases only by a factor of about 7. 



6 Conclusive remarks 

Let us list some key points of the results achieved: 

(1) The starting point in the analysis was the Gott-Hiscock exterior metric 
(5), given for the annular region a < r < R. All three types of the classical 
modes were worked out: The TM, TE, and TEM modes. The calculation 
was based upon perfect conductor boundary conditions at r = a, R. 

(2) Quantum mechanically, only the TM modes were considered, cor- 
responding to the restriction (12). The continuity conditions on the field 
component at the instant t — to were Eqs. (48) and (50). The energy 
production was found to be strongly peaked in the lowest (m, s) — (0, 1) 
mode; cf. Eq. (65). The much lower energy production in the next mode, 
the (1,1) mode, was given in Eq. (66). This expression shows the charac- 
teristic property of being proportional to {G/i/tY; this is the same kind of 
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behavior as found in the same mode in |^ and also for the scalar field in the 
azimuthally symmetric mode in p[. 

(3) The mathematical reason for the differences from the zero-width case 
lies in the need to include the Neumann function in the solutions for the 
mode functions; cf. Eqs. (43) and (15). The value of the energy produced 
was found to be greatly influenced by the roots of the boundary conditions, 
for t < to, and for t > to- 

(4) We assumed a GUT-scale string, and the input parameters (string 
mass, radius, waveguide dimensions), were chosen accordingly. The formal- 
ism of course permits scaling of these parameters. Inserting various values of 
the ratio a/Rwe found that the value of the produced energy did not vary 
very much (less than about one order of magnitude) when the variations took 
place in the range allowed by the approximations used. 



A On the evaluation of / 



ss' 



In the expression (58) for J^j,' we approximate g^^s by Pms, where the latter 
quantity is a root of J\m\{PmsR) = 0. This corresponds to the second term 
in the condition (12) for the TM mode being negligible. Moreover, setting 
f]m ~ m the integral in Eq. (58) becomes 



J\m\{Pmsr)C\m\{pms'r)rdr = - r^J\m\+i{pmsr)C\m\+i{Pmsr 



R 



(Al) 



where the orthogonality in s stems from Eq. 11.3.29 in |T^. Remembering 
that a/R and making use of N\m\+i{pmsR) / N\m\{pmso) — as well as 



iVi 



\m\+l 



(^) 



iVi 



\m\{z) 



(|m| + 1)- 



(A2) 



(cf. Eq. 9.1.9 in (1151), we see that the contribution from the lower limit 
r = a in Eq. (Al) can be omitted. 

Moreover, the second factor in the denominator in Eq. (58) can be written 

1 R 



r {J| 



J\m\ i^Pms^^ 



\r I \^ Iml + l 

>|m|(Pmsa) 

SO that the expression (58) is approximately equal to 



(A3) 



(A4) 
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as given in Eq. (59). As a corollary, we carried out numerical checks on the 
expression for Igs' using Maple V, supporting the result above. 
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(m,s) 






(0,1) 


2.40483 


2.54821 


(1,1) 


3.831711362 


3.831711363 


(1,2) 


7.015592430 


7.015592434 


(2,1) 


5.135632452 


5.135632452 



Table 1: Roots for various modes {a/R — 10 ^) 



a/I 


l 




Xf3U 


10- 


-1 


3.313938715 


3.940945907 


10" 


-2 


2.800921755 


3.832889531 


10- 


-3 


2.654814168 


3.831723171 


10- 


-4 


2.587120412 


3.831711480 


10- 


-5 


2.548210049 


3.831711363 



Table 2: Roots for various values of a/R 



a/R 


Woi X ALR 




Wn X ALR 




10-1 


2,493918307x10- 


-1 


3.026897060x10- 


-2 


10-2 


5.601176720x10" 


-2 


3.527890658x10- 


-7 


10-=^ 


2.353818082x10- 


-2 


3.229962982x10- 


-12 


10-^ 


1.284358731x10- 


-2 


6.036824072x10" 


-11 


10-5 


8.066977810x10- 


-3 


6.130061838x10" 


-11 



Table 3: Energy production in the modes (0,1) and (1,1) 
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